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,__! Abstract. A model for the transmission of dengue disease is presented. It consists of eight mutually-exclusive compartments 

representing the human and vector dynamics. It also includes a control parameter (insecticide) in order to fight the mosquitoes. 
The main goal of this work is to investigate the best way to apply the control in order to effectively reduce the number of 

qv infected humans and mosquitoes. A case study, using data of the outbreak that occurred in 2009 in Cape Verde, is presented. 
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O INTRODUCTION 

Dengue is a mosquito-borne infection mostly found in tropical and sub-tropical climates worldwide, mostly in urban 

Sand semi-urban areas. It can cause a severe flu-like illness, and sometimes a potentially lethal complication called 
dengue haemorrhagic fever. About 40% of the world's population are now at risk. The spread of dengue is attributed to 
expanding geographic distribution of the four dengue viruses and their mosquito vectors, the most important of which 
^_l is the predominantly urban species Aedes aegypti. The life cycle of a mosquito presents four distinct stages: egg, larva, 

^- pupa and adult. In the case of Aedes aegypti the first three stages take place in or near water while air is the medium for 

ON the adult stage. The adult stage of the mosquito, in the urban environment, is considered to last an average of eleven 

**l days. 

The paper is organized as follows. The next section presents a mathematical model of the interaction between human 
and mosquito populations. Then, the numerical experiments using different strategies for the insecticide administration 
are reported. In the last section the conclusions are presented. 
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THE MATHEMATICAL MODEL 



Considering the work of [4], the relationship between humans and mosquitoes are now rather complex, taking into 
account the model presented in [2]. The novelty in this paper is the presence of the control parameter related to adult 
mosquito insecticide [6], 

The notation used in our mathematical model includes four epidemiological states for humans (index h): 

Sh{t) susceptible (individuals who can contract the disease) 
Ef,(t) exposed (individuals who have been infected by the parasite 

but are not yet able to transmit to others) 
7/j(f ) infected (individuals capable of transmitting the disease to others) 
Rh(t) resistant (individuals who have acquired immunity) 

It is assumed that the total human population (N/,) is constant, so, N^ = 5/, +£/,+//, +Rh- There are also other 
four state variables related to the female mosquitoes (index m) (the male mosquitoes are not considered in this study 
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because they do not bite humans and consequently they do not influence the dynamics of the disease): 

A m (t) aquatic phase (that includes the egg, larva and pupa stages) 
S m (t) susceptible (mosquitoes that are able to contract the disease) 
E m (t) exposed (mosquitoes that are infected but are not yet able 

to transmit to humans) 
I m (t) infected (mosquitoes capable of transmitting the disease to humans) 

In order to analyze the effects of campaigns to combat the mosquito, there is also a control variable: 

c(t) level of insecticide campaigns 

Some assumptions are made in this model: the total human population (N/,) is constant, which means that we do 
not consider births and deaths; there is no immigration of infected individuals to the human population; the population 
is homogeneous; the coefficient of transmission of the disease is fixed and do not vary seasonally; both human and 
mosquitoes are assumed to be born susceptible; there is no natural protection; for the mosquito there is no resistant 
phase, due to its short lifetime. 

The Dengue epidemic can be modelled by the following nonlinear time-varying state equations for Human Popula- 
tion and vector population, respectively: 



d ^{t)= i l h N h -(Bp mh I ^ h 



V-h)Sh 



^(t) = n b (l-^)(S m + E, 



<%(?)= B&, 



Im < 



dt VJ — "h>mh J^Sh — ( Vh + Hh)Eh 

d k{t) = v h E h -{t] h +Hh)h 

*lhh — VhRh 



(t) = -(Bfa 



I 



■#« 



dt 

ds„ 

dt 

*fr{t)=Bp hn fc£ 

dim 

dt 



(0 



'Im^'tf 



' f^rn'm 



n+Im)-(riA-\-LlA)A n 
l^mj^m ~r ^}A^m Co m 

Clin 



(1) 



with the initial conditions: S h (0) = S h0 , E h (0) = E h0 , 4(0) = 4o, ^/,(0) = Rho, A m (0) =A m0 , S m (0) = S m0 , E,„(0) = 

E,n0, 4l(0) = InO- 

Notice that the equation related to the aquatic phase for the mosquito does not have the control variable c, because 
this kind of insecticide does not produce effects in this stage. 
The parameters used in the model are: 



Nfi total human population jj.a 

B average daily biting (per day) tj^ 

P m h transmission probability from I m (per bite) l/^m 

Phm transmission probability from 4 (per bite) 1/v/, 

1 l\ih average lifespan of humans (in days) m 

1/rjh mean viremic period (in days) k 

\/fX m average lifespan of adult mosquitoes (in days) K 

jj,/, number of eggs at each deposit per capita (per day) 



natural mortality of larvae (per day) 
maturation rate from larvae to adult (per day) 
extrinsic incubation period (in days) 
intrinsic incubation period (in days) 
female mosquitoes per human 
number of larvae per human 
maximal capacity of larvae 



The Figure 1 shows the relation between human and mosquito and the corresponding model parameters. 







FIGURE 1. Interaction between human and mosquito 



NUMERICAL EXPERIMENTS 

The values related to humans originate from data collected in an infected period in Cape Verde [1]. However, as it was 
the first outbreak that happened in the archipelago, it was not possible to collect any data for the mosquito. Thus, for 
the Aedes aegypti we have selected information from Brazil where dengue is already a long known reality [8, 9]. 

The numerical tests were carried out using Scilab [7] with the following values [8, 9]: Nk = 480000, B = 1, 
A,,,, = 0.375, &„, = 0.375, fl h = 1/(71 x 365), t] h = 1/3, Mm = 1/11, lib = 6, fi A = 1/4, T] A = 0.08, t] m = 1/11, 
V/, = 1/4, m = 6, k = 3, K = k X Nf,. The initial conditions for the problem were: S/,o = Nf, — £/,o — 4o, £/i0 = 216, 
7/iO = 434, Rho = 0, A„,o = k x A'/,, S m o = m x A^, E m o = 0, 7„,o = 0. The final time was tf = 84 days. 
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FIGURE 2. Infected mosquitoes and individuals 

In the literature it has been proved that a disease-free equilibria (DFE) is locally asymptotically stable, whenever a 
certain epidemiological threshold, known as the basic reproduction number, is less than one. The basic reproduction 
number of the disease represents the expected number of secondary cases produced in a completed susceptible 
population, by a typical infected individual during its entire period of infectiousness [3]. In a recent work [5] it was 
proved that if a constant minimum level of insecticide is applied (c = 0.084), it is possible to maintain the basic 
reproduction number below unity, guaranteing the DFE. 
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FIGURE 3. Infected mosquitoes and individuals vs periodicity 

To solve the system (1), in a first step, several strategies of control application were used: three different frequencies 
(weekly, bi-weekly and monthly) for control application, constant control (c = 0.084 in [5]) and no control (c = 0). 
The weekly/bi-weekly/monthly frequency means that during one day per week/bi-week/month the whole (100%) 



capacity of insecticide (c = 1) is used during all day. Constant control strategy (c = 0.084) consists in the application 
of 8.4% capacity of insecticide 24 hours per day all the time (84 days). In this work, the amount of insecticide is an 
adimensional value and must be considered in relative terms. 

Figure 2 shows the results of these strategies regarding infected mosquitoes and individuals. Without control, the 
number of infected mosquitoes and individuals increases expressively. 

Realizing the influence of the insecticide control, further tests were carried out to find the optimum periodicity of 
administration which, from gathered results, must rest between one and two weeks. The second phase of numerical 
tests, Figure 3, considers four situations: 7 days, 1 1 days, 12 days and continuously c = 0.084. To guarantee the DFE, 
the curves must remain below the one corresponding to c = 0.084. 

The amount of insecticide, and when to apply it, are important factors for outbreak control. Table 1 reports the total 
amount of insecticide used in each version during the 84 days. 
TABLE 1. Insecticide cost 

7 days 11 days 12 days 15 days 30 days c = 0.084 

insecticide amount 12 8 7 6 3 7.056 



CONCLUSIONS 

The numerical tests conclude that the best strategy for the infected reduction is the weekly administration, however it is 
also the most expensive one (insecticide cost). The best result obtained is between 1 1 and 12 days, with the insecticide 
amount in the closed interval from 7 to 8, confirming the amount for c = 0.084 in [5], The 11 or 12 days between 
applications can be directly related to the span of adult stage for the mosquitoes, an average of eleven days in an urban 
environment. 

In this work the insecticide administration time was considered continuous in each day (24 hours per day). As a 
future work, an optimization problem will be formulated and solved in order to find the best plan in terms of periodicity 
and duration (limited number of hours per day). The plan must consider the advantages to apply the insecticide only 
during the night. 
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